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I. INTRODUCTION 

Dileptons have been proposed in the late 70's as penetrating probes [0] of hot and 
dense nuclear systems. They are presumably created in all stages of heavy ion reactions 
by several distinct mechanisms. Once produced, they practically do not interact with 
the surrounding hadronic matter. Low mass dileptons are of particular interest (see e.g. 
|0-[^). They can reveal information of the hadronic properties in the reaction zone. 
Several experiments have focused on low mass lepton pairs: the DLS spectrometer at 
the BEVALAC @,|], the CERES and HELIOS detectors at the SPS at CERN [0. The 



dilepton spectrometers HADES at SIS at GSI and PHENIX at RHIC in BNL [|T2 
are under construction. The most striking result of the high energy dilepton programs so 
far is the observed enhancement in heavy systems at low invariant masses as compared 
to "conventional" hadronic cocktails and models. A dropping mass [Il3i|l4| or dissolving 



spectral function 1T^|T^ of the p meson have been offered in attempt to explain these 
data. Recently, a systematic measurement of dilepton production at BEVALAC energies 
has been published. Data are available for elementary pp and pd collisions as well as for 
nucleus- nucleus collisions [HQ. 

The aim of the present work is to investigate dilepton production within the micro- 
scopic n-body transport model UrQMD. The outline of the paper is as follows: In section 
in a brief survey of the UrQMD model is given. A more detailed description of the im- 
plemented dilepton production mechanisms follows. Section |ITI| contains the calculations 
of the elementary pp and pd dilepton cross sections in comparison to recent DLS mea- 
surements. Section |I^ shows mass spectra for heavier systems. Summary and concluding 
remarks are given in section |V|. 

II. PRODUCTION OF RESONANCES AND DILEPTONS 



A. Hadron production in the UrQMD model 



The UrQMD model is based on the quantum molecular djTiamics concept [[r7| -p!9 
The Hamilton's equations of motion are solved for gaussian wave packets. The model 
allows for the production of all established meson and baryon resonances up to about 
2 GeV with all corresponding isospin projections and antiparticle-states. The collision 
term describes particle production by resonant excitation channels and, for higher masses, 
within a string fragmentation scheme. The UrQMD model is designed to cover consistently 
the whole range of bombarding energies per nucleon from 200 MeV to 200 GeV. For 
dilepton production at BEVALAC energies, the resonant production of neutral mesons is 
most important. A detailed description of the model can be found in refs. [^|2l |. 



The formation of light mesons at low energies is modelled as a multistep process that 
proceeds via intermediate heavy baryon and meson resonances and their subsequent decay 
p^ . The resonance parameters (pole masses, widths and branching ratios) are taken from 



23[] , but large uncertainties of these parameters are used to obtain a consistent fit to cross 
section data. For example, the production of uj mesons is described in the UrQMD model 
by the formation and the decay of the A^*(1900) resonance. It decays to 35% into A^vr 
and to 55% into Nuj. As suggested in ref. |2^, the rj production proceeds not only 
via A^*(1535), but invokes also nucleon resonances with masses 1650, 1700, 1710 and 
2080 MeV. A full list of the UrQMD resonance parameters is published in ref. [pO|] . 

Broad resonances have mass dependent decay widths. For the resonance cross sections 
in baryon-baryon, meson-baryon and meson-meson collisions conventional Breit-Wigner 
parametrisations are used with mass dependent widths. At higher energies, the resonant 
particle production does no longer predict the observed cross sections. There, the string 
picture, described in [PT|, is employed. 



Fig. ^ shows calculations of exclusive {pp — ^ mpp) and inclusive {pp — > mX) cross sec- 
tions for the production of neutral mesons m = vr^, r], p°, u;. Recent data on the exclusive 
1] production just above threshold are reproduced reasonably well [ES. Note that above 



3.5 GeV the exclusive cross sections become less important, because the string channels 



open and allow for multiple resonance production. For the present article the relevant y/s 
values are below 4 GeV, where so far no data on the inclusive channels are available. 

B. Dilepton radiation in UrQMD 

Dileptons can be produced in hadronic decays and collisions. The mechanisms that 
are expected to dominate in the low mass region (with invariant dilepton masses below 
1 GeV) are the Dalitz decays {A -^ Be^e~) of neutral mesons 7r^,ri,ri',uj, as well as 
the Dalitz decay of the A(1232) resonance. Around their pole masses the direct decays 
{A -^ e^e~) of the vector mesons p^^uj and are expected to dominate the spectrum. 
These modes are of special interest: the invariant mass of the dilepton equals the vector 
meson mass in the hot and dense medium. Also included in the present model calculation 
is the pn bremsstrahlung. It contributes especially at low bombarding energies. In accord 
with ref . [Q , other sources like Dalitz decays of heavier resonances, pp, ttN bremsstrahlung 
etc. are neglected for the reactions considered here. Direct dilepton production in binary 
collisions, e.g. tt+tt^, i^+iT^ -^ e+e^ or np -^ e+e^, is not evaluated. In the framework 
of the present model this would correspond to a partial double counting: these reactions 
are supposedly included as explicit multi-step processes (e.g. np ^ cj) ^ e+e^). This 
also holds for the Dalitz decay of the ai ^ np ^ vre+e^. At higher bombarding energies, 
however, direct processes like Tip —>■ 7re+e^ might significantly contribute to the dilepton 
yield [^ . 

1. Dalitz decays of mesons 

Decays of the type A -^ Be^e~ are not real multi body decays (for which usually the 
name Dalitz is used) but can be reduced to a decay into B plus a virtual photon (with 
an invariant mass M) and subsequent conversion of the latter. Thus the matrix element 
factorises [p7| , p^ 



\M\' = \M{A -. 57*)P^I^(7* - e+e- 



This reflects in the differential decay rate, which can be written as a product of the 
conversion rate of the virtual gamma [^ and the decay width TA-tBr* of ^ i^^o a massive 
photon: 



dVA^Be+e- _ 2a Ami L ^ 2ml 



l-—^ 1 + T7f TA^Br- (2) 



Here me is the lepton mass. If one assumes that A decays isotropically in its restframe 
the A — > i?7* width is given by 

^A-^B-y* = Q 2"I-^A^-B7*| • (3) 



iTcm, 



Here |pcm| is the decay momentum 



X2(m\,m%,M'^) 



\Pcm\ = -^^^^ ^ (4) 

with the kinematic function X{x'^ , y'^ , z^) = (x^ — {y + zY){x^ ~ iv ~ ^)^)- 

A remaining difficulty resides in the calculation of the matrix element of the A —>■ Bj* 
decay. The meson decays considered here are either of the type a) pseudoscalar (vr", r], r]') 
into a vector particle (7) plus a virtual photon or b) vector meson (u) into a pseudoscalar 
meson plus a virtual photon. In both cases one gets ||3^ 



\MA^Br\' = l\fAB{M^)\'\{ml,ml,M^) . (5) 

The form factor JabIM^) is introduced to account for the strong interaction part of the 
vertex. It is common to normalise to the decay width into real photons [BO] by dividing 



{mA — m^) 2 /«\ 

r^^B. = ^^^^^^^l/AB(o)| , (6) 



which re-expresses the form factors to FABiM"^) = f(Q)- 

The form factors can be obtained from the vector meson dominance model (VMD). 
In the present calculations the following parametrisations are employed p8|,|6| 



F^o(M2) = 1 + b^oM^ 



F.(M')-[i a; 



A/'\"' 



F^M^ 



F^i{M'^ 



(A2_M2)2 + A2^2 



(A2, _ M2)2 + A2,^2, 

with 6^0 = 5.5GeV~^ A^ = 0.72 GeV, A^ = 0.65 GeV, 7^ = 0.04 GeV, A^^/ 
and7; = 0.10GeV. 



(7) 
0.76 GeV 



2. Delta Dalitz decay 

The situation is more comphcated for the Dahtz decay of the A(1232) resonance. 
To complete eqs. (^ and (|^), the matrix element for the process A -^ N'-f* has to be 
calculated. The corresponding interaction vertex has been analysed by Jones and Scadron 



| 3T[| in the form 



Cint = ef^«r^M.^ + h.c. 



(8) 



where ip, \l//3 and A^ represent the fields of the nucleon, of the delta and of the photon 
respectively. The dominant magnetic dipole transition yields the vertex function 

-3(mA + mAr) / fl„ n,, 1 



r^M = Gm{M^ 



Ba . 0a . ^ 0a 

-mAXi 75 + X2 75 + 7:X3 75 



(9) 



2m^((mA + mjv)2-M2) V "'^'"" " ' "^ " ' 2 
The choice of the x bears some freedom if only current conservation is ensured (g/jF^^ = 0). 



Following ref. ||31|, one may write 



xf 



{q^r -q-19^ 
(q^P^" -q-Pg 



q^q^" - q^g^^ 



(10) 



with P = 1/2(pa + Pn) and q = Pa — Pn- Gm in eq. @ stands for the magnetic 
dipole form factor. It can be fixed at M = to the decay into a real photon, yielding 
Ca/ (0) = 3.0. The M-dependence of the overall form factor is subject to speculation: the 
time-like electromagnetic form factors for baryons are unknown in the kinematic region of 



interest. In the present calculations any VMD-type form factors are omitted. This gives 
a lower limit for the A Dalitz contribution in the region of the vector meson poles. 
Using eqs. (|,^, we can express the matrix element via 

Mi = eUp{pA, SA)Xi^75(^ti{P-y, S^)u{pn, Sn) ■ (H) 

The matrix element may be written as a linear combination of the Aii and Ai*: 

where ci = —rriA, C2 = 1 and C3 = 1/2. The appearing traces have been calculated using 
the Mathematica package HIP [^. One then obtains: 
3 



^ « 1 
^ CiMiCjM* = - ((mA - rriNf - M^ 



i,j='^ 



X 



(jml + Um\M'^ + 3M^ + 8m^m^ + 2m\m% + 6M^m% + 3m^) . (13) 



Substituting eqn. (|T2D and (|1^) into eq. (^, we get the A Dalitz decay width. 

Fig. shows the differential mass distributions of the Dalitz decay probabilities as im- 
plemented into the UrQMD model. While the very low masses are dominated by 7r° and 
rj decays, the 00 and rj' decays are more important at higher masses. The contributions 
of the A decays are also shown for different masses of the A(1232). Heavy A's natu- 
rally contribute more to the probabilities than lighter ones. In calculating the dilepton 
spectrum these probabilities are multiplied with the yields of the corresponding particle 
species. At low energies the dominant vrA^A system [|18[ can push the A contributions 



over those of the heavy mesons. On the other hand, the small yield of the rj' will make it 
invisible, due to the tremendous background of other sources in the relevant mass region. 

3. Direct decays of vector mesons 

The decay width of the electromagnetic two-body decays of the vector mesons are 
assumed (similar to ref. Q) to be of the form 



V-...-(M) = iij^l - M A + ?^J e(M - 2m.) . (14) 



The constants Cy are fitted to yield the vacuum widths given in ref. [^ at the resonance 
poles. One obtains Cp = 3.079 x 10"^ GeV^ c^ = 0.287 x lO^^ GeV^ and c^ = 1.450 x 
10-*^ GeV^ 

It is assumed that vector mesons can radiate off dileptons continuously [Q. The 
dilepton mass distribution is given by the time-integral over the mass distributions of 
mesons plus an additional term which accounts for the decays after some typical freeze- 
out time tf (the last term can be dropped if tj — > oo) 



dM Jo dM " """^ ' dU ^v,tot{tf) 

If absorption is negligible, this approach is equivalent to the method of adding one dilepton 
(with appropriate normalisation) at each decay vertex. However, this "shining" method 
gives a better sampling of the density probed by the hadron and thus a better statistics for 
density dependent spectral functions. In addition, there is generally a complicated time 
dependence of the total hadron width rtot(t,r), caused by the in-medium modifications 
of the quasiparticle widths. The approach eq. (|15D has the advantage to be explicitely 
independent of the in-medium width, if tj is reasonably high, so that the total width takes 
on it's vacuum properties. 

^. Incoherent bremsstrahlung 

In the soft photon approximation (SPA) the cross section for real photons with four- 
momentum g'^ in collisions a + b -^ X can be expressed as 

dcr^ denotes the differential cross section of the (strong) ab interaction without any photon 
in the final state, e is the polarisation vector of the photon. The current J^ is given by 

Pa-q Pb-q ~i Pi-q 



where the Q's and p's denote the charges and momenta of the corresponding particles. 

The SPA is justified only for M, go —* 0. To extrapolate to the case of hard and 
massive virtual photons, a phase space correction can be applied by multiplying the cross 
section with the ratio of the phase space integrals with/without a virtual photon [Q. 
Similarly to eq. (0) one gets 



(PqdM eTT^V M2 1^ M2y' ' Rnis)qoM' 

Here R„ is defined as 



i-n 



Rn{s) = J d^nis,Pi...Pn) , (19) 

where d^n is the volume element of the n-dimensional Lorentz invariant phase space and 
s is the squared effective energy of the system after the emission of the 7*, 

s = s + M^ - 2v^go • (20) 



The correction factor for two outgoing particles reads 



(21) 



i?2(s, ml, ml) Ai/2(s, ml, ml) s ' 

Eq. ([TsD is a general expression for the bremsstrahlung dilepton production in the 
SPA. In the case of proton-neutron bremsstrahlung {pn — > p'n''~f*), the current is given 
by J^ = Pp//q ■ Pp' — Pp/q ■ Pp- As a result we have 

\^-jf = ^, (22) 

q^mj 

with the Mandelstam variable t = (pp — ppi)'^ . Then eq. ([I8| ) takes the form 



rf<r"^- _ a' J7~^ (. I 2mn cT R,{-s) 
d^qdM Q7i^\ M2 1^ M^)qlMR2{s)' ^ ' 

Here a is the momentum-transfer weighted pn elastic cross section 

a = r (4) pt . (24) 

J-(s-4m2) \m^J dt 



A proper parametrisation of a is rather important [^. A systematic overestimation of 
the differential bremsstrahlung production results if the asymmetry of the momentum- 
transfer weighted cross section is not included. In the present work, a parametrisation 
similar to that of ref. [^ is used. It consists of a symmetric low energy part and an 
asymmetric high energy component. 

Fig. D shows the cross section for pn —^ pn + e~^e~ bremsstrahlung at different bom- 
barding energies. Note that at low masses, the cross section is practically not sensitive 
to the incident energy. As a result, the pn bremsstrahlung is relatively unimportant as 
resonance channels become dominant at E^m ~ 1.5 GeV. 

Bremsstrahlung from pp collisions is expected to be small at low energies due to 
destructive interferences. However, at the highest energies considered here, the pp 
bremsstrahlung contribution may be comparable to the pn bremsstrahlung yield [^ 



On the other hand, the total contribution of bremsstrahlung is found to be negligible for 
all practical purposes at those energies (see sect. |111BD . Note that the SPA already gives 



an upper estimate of the expected dilepton yield |37 



III. DILEPTON CROSS SECTIONS IN ELEMENTARY SYSTEMS 

Before one can investigate the dilepton production in nucleus-nucleus reactions, one 
should first check the model in elementary pp and pd collisions. Recently, the DLS col- 
laboration published a systematic study on dielectron cross sections in light systems for 
beam kinetic energies Ekin from 1 to 5 GeV. The corresponding centre-of-mass energies 
range from below the t] production threshold up to 3.6 GeV, where phase space is wide 
open for the abundant production of a large variety of resonances. 

The acceptance of the DLS does not cover the entire phase space. Thus the calcu- 
lated dilepton pairs are corrected for the limited spectrometer acceptance region (DLS- 
filter V4.1). The finite mass resolution of AM/M ^ 0.1 is incorporated by folding the 
calculated spectra with a Gaussian of width AM. The acceptance correction strongly 
influences the low-mass spectra while the mass resolution smoothing affects the shape of 
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the spectra at higher masses. 

A. pp 

The resulting dilepton mass spectra for pp colhsions are shown in fig. ^. 

Only 0.46 GeV of cm. energy are available for particle production at Ekin = 1-04 GeV. 
Thus, only pion- and A-Dalitz decays contribute in this case. The model satifactorily 
reproduces the data at low masses, but underestimates them around Me+e- = 0.4 GeV. 
The disagreement could be caused by the abovementioned uncertainties in calculating the 
electromagnetic form factor of the AN'j system. The first generation of DLS data was 
incompatible to free form factors [^p8|. This situation is now unclear, because the data 
of the second run tends to exceed the first. However, limited statistics and the better 
agreement at higher beam energies precludes a definite conclusion on the origin of the 
"enhancement" at Ekin = 1.04 GeV. 

At Ekin = 1-27 — 1.85 GeV the model explains the data with growing infiuence of 
meson decays. One can see from fig. ^ that first the r^-Dalitz decay and then the direct 
p and u decay become more important. In the UrQMD calculation oi pp, most of the p^ 
(~ 99%) are not created in tt+tt^ collisions, but result from heavy baryon decays. The 
limited phase space at low cm. energies [^ and the strong (M~^) mass dependence of the 
dilepton decay widths (0) are responsible for the deformation of the p spectrum towards 
low masses. The cross section for the Dalitz production through A decays remains rather 
constant at low masses, with increasing energy, but increases towards higher masses. 
Indeed, it is known from pp -^ nA'^^ reactions, that heavier A's become more important 
with increasing beam energy due to the mass dependence of the decay widths |^. In 



fact, there is a good agreement at intermediate beam energies (disregarding the data point 
at Me+e- = 200 MeV) which would be worsened if the electromagnetic AA^7 form factor 
would be included, in particular at E^in > 1.6 GeV. 

The model does not reproduce the shape of the measured spectrum at E^in = 2.1 GeV 
and 4.9 GeV. Around M = 0.6 GeV the ratio between the pp and pd data decreases for 
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the two highest beam energies. If this gets confirmed, it might be explainable by some 
strong contribution of pp bremsstrahlung, which was not considered here. 

B. pd 

pd calculations are shown in fig. ^. The calculated events are minimum-bias triggered 
at a maximum impact parameter of 6 = 1.6 fm, i.e. a geometric cross section equal to the 
measured one of about 80 mb. 

An important difference to the pp system is due to the internal motion of bound 
nucleons. This motion smears out the production thresholds. Consequently, one observes 
subthreshold contributions from rj and p mesons already at Ekin = 1.04 GeV. 

Unlike in pp collisions, the proton-deuteron data overestimate the model results in 
the mass region dominated by the t] decay. This indicates an asymmetry in the pp and 



pn production cross sections as predicted from one-boson exchange models |^. Indeed, 
it was measured that just above the threshold, the ratio of pd to pp induced r] cross 
sections is much larger than two |^ . However, due to the Fermi motion in the deuteron. 



it is not straight forward to extract pn cross sections from pp and pd measurements. On 
the other hand, direct measurements of the rj production in pn are still not available. 
We have therefore evaluated the ratio i?^ = (7{pd -^ r]X)/a{pp —>■ rjX) without any 
explicit modification of the pn cross section to estimate the infiuence of Fermi motion. 
At Ekin = 1.27 GeV, close to threshold, we get the large value i?^ = 17 ± 3 while above 
Ekin = 1-61 GeV, Rff remains practically constant, between 2.8 and 2.3. 

However, the large value of i?^ is not enough to explain the integrated total dilepton 
ratios. This is shown in fig. ^ where 

R = ,f ^"^^.'^^^ (25) 

fMmax I do'pp \ ^ ]\/f 

J0.15 GeV \dM ) "^^-' 

is plotted vs. the incident energy. The predicted dilepton ratios are rather low at energies 
smaller than 2 GeV. A similar behaviour is also found in BUU calculations ||^, which 
include pp bremsstrahlung. This indicates that there is a common feature of known 
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transport models with "conventional" sources to underestimate the recent DLS nucleus- 
data at intermediate dilepton masses [§. At Ekin = 1.27 GeV the model predicts only 
R = 5.2. This value is lower than the corresponding one of i?^, because the rj is yet only 
a marginal dilepton source at this energy. 

Within the statistical uncertainties the mass-differential dilepton cross sections for 
pd can be explained by the model for the two highest beam energies. Note that pn 
bremsstrahlung is relatively unimportant. 

IV. DILEPTONS IN HEAVY ION COLLISIONS 

Fig. 1^ compares the UrQMD predictions with DLS data for various nucleus-nucleus 
reactions [§]. At first glance the examined systems all exhibit a qualitatively similar 
behaviour in three distinct mass regions: 

In the lowest mass region (up to 150 MeV) the spectrum is dominated by the pionic 
Dalitz decay. There are also strong effects of the acceptance filter, which significantly 
suppresses the low- mass yield (compare Fig. ^jfor the shape of uncorrected Dalitz spectra). 
Other sources are of little importance in this mass region. However, near M = 150 MeV, 
the A Dalitz decay becomes more important. 

At intermediate masses the data show a noticeable enhancement of the dilepton yield 
as compared to the model calculations. There is a considerable confusion about these 



recent data because the new yields strongly exceed earlier published measurements [^ 
The latter have been revised by the DLS collaboration due to large, previously uncorrected 
trigger inefficiencies. However, the present calculations and, also, results of other models 
e.g. P,^3|-|i5[], are closer to the earlier measurements in the considered mass range (see 

The theoretical spectrum at M > 600 MeV is dominated by direct p° and uj decays. 
The model cross sections for dilepton production via p mesons nicely account for the data 
in this region. The two highest data points in He-Ca suffer from lack of statistics. For all 
nuclear systems, only about 50% of all p mesons stem from tt+tt^ annihilations, the other 
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50% are produced in decays ||22[ . 



Turning back to the intermediate region, one sees a similar dilepton enhancement as 
in the pd data. There one may partially attribute the enhancement to the high pn -^ rjX 
cross section. Therefore, it seems reasonable to reproduce the data by artifically enhancing 
the 7] yield. We rescaled the rj yield by a factor /^ and found good agreement both in 
shape and absolute yield with /^ k, 10 for C+C and k, 20 for Ca+Ca. However, this is 
not supported by TAPS data which found the 47r extrapolated r] cross section in Ca+Ca 



at 1 GeV to be about 20 mb [^. Our rescaled cross section corresponds to ^120 mb. A 
dropping rj mass has been introduced in ref. |^ to explain the DLS data. But these studies 
also found the increased rj cross sections to be incompatible to the TAPS measurements. 

We note that the inclusion of a density-dependent AA^7 form factor for the A Dalitz 
decay, as discussed in ref. , can give a sizeable enhancement of the calculated diletpon 
yield at intermediate masses 0.2 GeV < M < 0.6 GeV. However, this is ruled out as an 
explanation, because the fair agreement of the calculation to the data in the high-mass 
region M > 600 MeV would be destroyed by a strong overestimation. 

Due to a — uj mixing via an NN~^ loop , the 2tt decay channel of the u might be 
significantly enhanced in nuclear matter |^^. To get an estimate for possible effects on 
intermediate mass lepton pairs, we have increased r^_,7r7r by a factor of 500 to approx- 
imately 100 MeV. The solid curves in fig. |^ show the result of this calculation for the 
Ca+Ca system. As can be seen, the mass distribution of dilepton pairs from direct u 
decays becomes very broad. Nevertheless, this contribution has almost no effect on the 
total dilepton yield and does not suffice to explain the measured enhancement. 

One approach to describe the density dependence of vector meson masses is the linear 
scaling law of Hatsuda and Lee ||48|| : 



m;^ = my(l-0.18^/^o) , (26) 

where f)(o) is the local (ground state) density and my is the (modified) vector meson 
mass. Although more sophisticated calculations predict a more complex behaviour of 



the vector meson spectral functions [15], the scaling law seems to be reasonable for the 
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effective masses. To ctieck tfiis idea using ttie UrQMD model, tfie poles of tlie produced 
quasi particle vector mesons have been shifted according to (PB]). One obtains the dashed 
curves in fig. p. It is found that the dropping mass scenario - and also more complex 
in-medium spectral functions (see ref. [^) - cannot account for the new DLS data. 



V. SUMMARY 

Dilepton production has been studied within the microscopic non-equilibrium trans- 
port model UrQMD. The production mechanisms have been critically revisited. We have 
compared the model with the DLS data for pp and pd collisions at different beam energies. 
Resonance decays into dileptons were found to be able to explain the low energy pp data. 

The UrQMD model predictions for the pd system are below the new DLS data for 
masses 0.3-0.6 GeV. A similar, but even stronger underestimation takes place in collisions 
of heavier nuclei. The present paper points out that the large enhancement in the data as 
compared to model calculations cannot be accounted for by two distinct hypothesises on 
the in-medium modifications of vector mesons. An enhancement of the rj production in 
pn collisions is able to reproduce the yield and shape of the AA spectra. However, huge 
Tj production cross sections would be required, in contrast to TAPS data. 
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FIGURE CAPTIONS 



FIG. |T] Cross sections for neutral meson production in pp collisions. Calculations are 
shown for the exclusive and inclusive production of 7r°, rj, uj and p° mesons in 
comparison to available data [^,^ • 



FIG. Differential probability distributions for Dalitz decays. The mass-differential 
branching ratios for decays of 7r° (dashed), rj (solid), uj (long-dashed) and rj' (dotted) 
mesons can be seen. Also shown is the parametrisation for Dalitz decays of the 
A(1232) resonances at pole mass (dotted) and of a heavy A(1232) with mass rriA = 
1.432 GeV (dash-dotted). 

FIG. Differential cross section for the production of a bremsstrahlungs pair in elastic 
pn collisions for beam kinetic energies E^in =1.04 (solid), 2.09 (dotted) and 4.88 GeV 
(dash-dotted). 
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FIG. ^ The e^e~ mass spectra for pp reactions at six kinetic energies. The upper sohd 
curve is the sum of all contributions. One can see the contributions of 7r° (dotted), 
A (dash-dotted), rj (dashed) and oj (dotted) Dalitz decays as well as of the direct 
p° (solid) and uj (thin-solid) dilepton decays. Every single dilepton has passed the 
DLS Filter 4.1 and a mass resolution of 10% is adopted. The data are from 0. 

FIG. ^ Acceptance-corrected mass spectra for proton -|- deuterium. The individual 
sources are marked are as in fig. ^, but there is an additional contribution from pn 
bremsstrahlung (thick-dotted). 

FIG. ^ Ratio of the integrated cross section for pd to pp reactions as a function of 
the beam energy. The solid line shows the UrQMD result with the DLS filter and 
resolution. The experimental data are taken from ref. 0. Only the larger systematic 
error bars are displayed. 

FIG. Model calculations for dilepton spectra from nucleus-nucleus collisions at beam 
energies of about 1 GeV in comparison to DLS data [Q. The systems d-l-Ca, a+Ca 
and Ca-|-Ca as well as C-f-C have been examined. See figs. |^ and ^ for additional 
information. 

FIG. ^ Total dilepton spectra from Ca+Ca collisions compared to two generations of 
data. To compare to the earlier measurements [Q the old DLS filter VI. 5 (dashed 
line) has been applied. The solid line corresponds to the calculation of fig. |^. 

FIG. Total dilepton spectra from Ca-|-Ca collisions with an increased uo width. Also 
shown are effects due to density dependent modifications of the vector meson masses. 
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Fig. I 
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Fig. g 
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